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ABSTRACT 

In this paper we present the Fermi All-sky Variability Analysis (FAVA), a tool to 
systematically study the variability of the gamma-ray sky measured by the Large Area 
Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. For each direction on 
the sky, FAVA compares the number of gamma rays observed in a given time window 
to the number of gamma rays expected for the average emission detected from that 
direction. This method is used in weekly time intervals to derive a list of 215 flaring 
gamma-ray sources. We proceed to discuss the 27 sources found at Galactic latitudes 
smaller than 10° and show that, despite their low latitudes, most of them are likely of 
extragalactic origin. 

Subject headings: keywords 



1. Introduction 

In 1844 the astronomer F. W. A. Argelander performed one of the first systematic studies 
of the variability of the night sky. He laid the study of variable sources "most pressingly on the 
heart of all lovers of the starry heavens, to perform an important part towards the increase of 
human knowledge, and help to investigate the eternal laws which announce in endless distance 



the Almighty power and wisdom of the Creator" ( Percy 2007 ) . Nowadays astronomers are not as 



poetic, but time has provided us with exceptional instruments for the quest. 

In this paper we present a systematic study of the temporal variations of the gamma-ray 
sky measured by the LAT on board the Fermi Satellite. The gamma-ray sky above 100 MeV is 
dominated by the Galactic diffuse emission, which originates from cosmic-ray interactions with 



interstellar matter and photon fields (e.g. Ackermann et al. 2012a). Additionally, an isotropic 



diffuse gamma-ray emission is detected, and it is the strongest source of emission at high Galactic 
latitudes ( Abdo et al.||2010e ). Both diffuse components are expected to be stable over the duration 



of the Fermi mission. On-top of this background, 1873 gamma-ray sources have been detected 
during the first two years of Fermi mission and reported in the second Fermi-LAT catalog (2FGL, 
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Nolan et al.|2012 ). Out of these sources ~24% are found to be variable on monthly time scales. The 
vast majority of the variable sources are associated with Active Galactic Nuclei (AGN), which are 
known to be variable across the electromagnetic spectrum (e.g. Aharonian et al.j|2009 Ackermann 
eTaLl|2011a|b[ ). 



Variability in gamma rays has so far been established only for a few sources in our Galaxy. 
Orbital modulation and isolated flares have been reported from 7 X-ray binaries, in which a neutron 



star or a black hole orbits a massive companion (Aharonian et al. 2005, 2006 Albert et al. 2006 



Abdo et al. 2009b; Hinton et al. 2009 Sabatini et al. 2010 Ackermann et al. 2012b). Variable 



gamma-ray emission has also been reported from the direction of Eta Car, a massive star-Wolf- 
Rayet binary (iTavani et al. 1120091 lAbdo et al.||2010cl iReitberger et all 120121). Recently, flaring 



gamma-ray emission has been found for two new source classes: Nova explosions (Abdo et al 



2010d) and the Crab pulsar wind nebula (Abdo et al. 2011c; Tavani et al. 2011 Buehler et al 



2012). The latter was thought to be a stable gamma-ray emitter, but was discovered to be flaring 



with the method described in this paper. In addition several other gamma-ray transients have been 



detected near the Galactic equator, but they are likely associated with distant AGN ( Vandenbroucke 
et al.||2010[ |Cheung et aT1|2012a ). 



The reasons for the detection of so few variable gamma-ray sources within our Galaxy remain 
unclear, whether due to astrophysical reasons, due to the statistically limited flux sensitivity, or 
due to systematic difficulties of detecting them owing to uncertainties in the modeling of the strong 
foreground of the Galactic diffuse emission. In this paper we present a new method developed to 
search for transients in the gamma-ray sky that does not require a diffuse emission model. We first 
describe the methods and proceed to assemble a list of flaring gamma-ray sources seen over the 
sky during the first 47 months of the Fermi mission. We then focus on the sources detected at low 
Galactic latitudes, as they may be of Galactic origin. 



2. The Fermi All-sky Variability Analysis 

Flux variability of LAT sources is usually studied with a maximum likelihood analysis, in 
which parameters of a model describing the point sources and diffuse gamma-ray emission in a 
given region of the sky are jointly optimized. The sensitivity of this approach is often limited 



by the uncertainties of the diffuse emission modeling, particularly in the Galactic plane (Nolan 



et al. 2012). Small inaccuracies in the instrument response functions can lead to time-dependent 



residuals which depend on varying observation conditions, e.g. off-axis angle of the sources or orbital 
position of the LAT, limiting variability studies ( Ackermann et aL|2012c ), A further limiting factor 
of the likelihood approach is that it is computationally intensive; it is currently difficult to perform 
variability studies in different time and energy windows over the entire sky. We therefore developed 
the Fermi All-sky Variability Analysis (FAVA), in which we search over a grid of regions on the 
sky for deviations from the expected flux based on the long-term average. While this approach is 
less sensitive than a likelihood analysis, it has three main advantages: 



1. The analysis is independent of any model for the diffuse gamma-ray emission. The diffuse 
emission is expected to be constant over the time of the Fermi mission. It therefore cancels 
out in the comparisons between number of expected and measured gamma-ray events. 

2. The analysis is computationally inexpensive, allowing us to blindly search for flux variations 
over the entire sky. The analysis is therefore unbiased, treating every direction on the sky 
equally, potentially yielding unexpected discoveries. 

3. No assumptions are made about the spectral shapes of the gamma-ray sources. Negative flux 
variations are treated the same way as positive ones. (Throughout this paper we refer to both 
positive and negative variations from the mean as flares.) 

We applied FAVA to the first 47 months of Fermi observations (2008-08-04 to 2012-07-16 
UTC), in weekly time intervals. The total number of weeks is 206. We considered two ranges of 
gamma-ray energy, E > 100 MeV and E > 800 MeV, to increase the sensitivity for spectrally soft 
and hard flares, respectively. We used the P7SOURCE_V6 event selection and only considered 
events with a zenith angle smaller than 95°, to limit contamination from the gamma-ray emission 
of the Earth Atmosphere, which is time variable in sky coordinates. 

We generate measured and expected counts maps with a resolution of 0.25 square degrees per 
pixel. The maps are smoothed by assigning to each pixel all events that were detected within a 
distance corresponding to the 68% containment radius of the Point Spread Function (PSF). The 
pixel positions are characterized in spherical coordinates by <j) and 8. The number N exp ((fi,6) of 
expected events in one pixel is derived from the number N tot ((j), 9) of events observed from the same 
direction over the first 47 months of observations. As the PSF depends on the photon energy E 
and on the incidence angle a with respect to the LAT ( Ackermann et aT]|2012c ), we integrate over 
these parameters: 



V^ V^ ATtotfi n\ ., i,j 



eff k ((f>,6) 



N^(<p,9)= ^ ^ N%{M)x gj— , (1) 

E:j=l..l2a:i=l.A *>J 

where e w and e tot are weekly and total exposures. We proceed to calculate the probability that the 
observed counts are a statistical fluctuation of the expected value, based on Poisson statistics. These 
probabilities are then converted to significance in units of Gaussian sigmas for easier visualization. 
The procedure is illustrated in Figure [T] for one example week in 2009 February, during which the 
first flare from the Crab Nebula was seen in LAT data. 

The sensitivity of FAVA to detect flares varies with the position in the sky, due to the anisotropy 
of the diffuse emission backgrounds. Additionally, the sensitivity depends on the energy spectrum 
of the flaring gamma-ray source. Typically gamma-ray sources have photon indices between 1.5 
and 2.5 in the Fermi-LAT energy range (Nolan et al.|2012>. The sensitivity for both of these cases 



are shown in Figure [2] Flares with a photon index greater than ~2 are typically detected at higher 
significance in the low-energy maps, while those with a smaller index are detected more significantly 
in the high-energy maps. 



The accuracy of FAVA was tested on simulations of a constant sky for 36 months of obser- 
vations. The simulated sky was composed of the Galactic (gal_2yearp7v6_v0.fits) and isotropic 
(iso_p7v6source.txt) diffuse emission modelaH and point sources. The latter were generated with 
random coordinates on the sky and according to the flux distribution derived in|Abdo et al. (201 Of). 



The distribution of the significance a of flux variations is displayed in Figure [3] for the simulations 
and the real data. In the low-energy interval, on average 154 events are recorded per sky pixel in 
each week. The simulated significance distribution is therefore expected to be close to Gaussian in 
the low-energy band. Indeed, a fit to the simulations shows that a Gaussian model with a standard 
deviation of 1.00 and a mean of 0.06 describes the distribution well, as shown in Figure [3J In the 
high-energy band an average of 3.6 events are recorded per sky pixel. The simulated distribution of 
flux variations is therefore not expected to be Gaussian. In particular, typically only a few counts 
are detected in several pixels at higher Galactic latitudes, leading to small-scale structure around 
|o"| < 1 in the significance distribution. A disagreement is visible in this region between the data 
and simulations. However, for |<r| > 2 the qualitative agreement with expectations is good also in 
the high-energy band, we therefore expect no significant biases in the detection of flares at high 
significance. As at low energies, one can see that the difference between simulated and real signifi- 
cance distributions increases with increasing a, corresponding to the real flux variations emerging 
above the statistical background. 

The good agreement of the statistical fluctuations in the low-energy band of the simulated 
constant sky with expectations shows that systematic effects are small in the simulations. However, 
as seen for the discrepancy in the |<r| < 1 interval in the high-energy band, additional effects might 
be present in the real data, as the simulations do not take into account possible sources of systematic 
errors. In particular, they do not account for any disagreement between the simulated and real 
instrument responses, nor for varying background levels due to residual cosmic rays mistakenly 
classified as gamma rays. It is difficult to assess these systematic effects from the data, as e.g. 
small flux variations might be present throughout the gamma-ray sky due to variable background 
sources. However, we can set upper limits on possible systematic errors by looking at presumably 
constant sources such as pulsars. Analysis of the brightest ones, the Vela and Geminga pulsars, 
shows that their relative count variations are compatible with a steady flux within <5% on weekly 
time scales. Systematic errors of FAVA for relative flux variations of bright sources are therefore 



comparable to those for the standard Fermi-LAT analysis (Ackermann et al. 2012c). 



3. List of flaring gamma-ray sources 

After calculating the significance maps for all weeks, we scan them for significant flares. To 



determine their positions we use the peak finding algorithm described in Morhac et al. (2000). 



Based on these detections, we build a list of flaring sources. For this we only consider flares with 



1 http://fermi.gsfc. nasa.gov/ssc/data/access/lat/BackgroundModels. html. 



significances greater than 5.5 a in the low-energy or high-energy band. This threshold was chosen 
so that the number of false flare detections due to statistical fluctuations is expected to be ~ 1 over 
the 206 weeks that were analyzed^ No flares were detected in the simulations of three years of data 
for a steady sky above this threshold. The number of false flares expected for 47 months of data is 
therefore <3 at 90% confidence. Additionally, we only consider flares that occurred far away from 
the average position of the Sun in the corresponding week. The Sun is a bright gamma-ray source 



that moves along the ecliptic by ~ 7° per week (Abdo et al. 2011b). We therefore only considered 



flares at a distance from the Sun >12° and >8° in the low- and high-energy bands, respectively. 
Finally, we merge low- and high-energy flares detected in the same week, if they are coincident in 
position within 3°, relating them to the same flaring event. For the position of the latter we use the 
position of the high-energy detection, due to its higher accuracy, as will be discussed in the next 
paragraph. A total of 1419 flares that fulfill the mentioned criteria were detected. Out of these, 
645 and 175 are detected in the low- and high-energy interval only, respectively. The remaining 
599 flares are detected simultaneously in both energy bands. 

To estimate the position accuracy achieved by the peak finding algorithm, we analyzed the 
distribution of flares around known flaring gamma-ray sources. As a source sample we chose the 
249 sources flagged as most variable in the second Fermi-L AT catalog (a variability index >83.2 
Nolan et al. 2012). The resulting distribution of flares per solid angle is shown in Figure Hi We 



m 



assumed that the reconstructed position of the peak finder follows a Gaussian distribution plus a 
constant background term from flares not associated with the sources in this representation. The 
best-fit model shown in Figure [4] represents the data in good approximation. We proceeded to 
calculate the distances within which 68% of the flares are contained in this parametrization. At 
low energies the radius is 0.8° and at high energies it is 0.6°. We verified that no systematic offset 
is present in the position estimation in any coordinate direction. 

Most variable gamma-ray sources, such as AGN or X-ray binaries, are known to have recurring 
flares. We therefore group the detected flares, and associate closely located flares to a single common 



flaring source. For this purpose we used a Minimum Spanning Tree (MST, |Nea et alT 2001). We 



first group the flares detected at high energies because their positions are better determined. We 
build the MST for these flares and merge neighboring flares with a distance of less than 2° in the 
spanning tree. We proceed to associate flares detected only at low energies to the ones found at 
high energy if their distance is less than 3°. Finally, we build the MST of the low-energy flares that 
were not associated and merge neighboring flares with a distance of less than 3° in the MST. 

The position of a FAVA source is found by averaging over the positions of its flares. If the 



2 The number of trials can be approximately estimated as the total sky area divided by the area of the PSF. Above 
100 MeV the 68% containment radius of the event-averaged PSF of the Fermi-LAT is ~ 3°. We therefore have 
~ 41253/ (it x 3 2 ) x 206 = 300558 independent tests in the sky. This results in ~0.01 expected false positives above 
5.5cr. For the high-energy maps the average PSF is ~ 0.6°, resulting in ~0.3 expected false positives above the same 
threshold. 



source is detected at high energies we use only the positions of these flares due to their better 
positional accuracy; otherwise the low-energy flares are used. The position error of the source is 
obtained by propagating the positional errors of the flares included. In addition to this statistical 
error, there is a systematic error, that can arise from false associations of flares to a source, as 
well as the finite binning of the sky maps. We estimated this error to be smaller than r sys = 0.1° 
by comparing the position of the FAVA sources to those in the list of variable 2FGL sources used 
previously. We assume r sys to be the systematic error on the source positions. 

A total of 215 sources are detected by FAVA. Out of these, 33 are detected at low energies 
only. Flares related to negative flux variations from the average emission are found for 22 sources, 
often during periods of quiescent emission. All of the latter also showed positive flares. No source 
was found which flared only due to a negative flux variation. Each FAVA source is referred to by 
its identification number composed of the right ascension in hours and minutes and the declination 
in degrees of the source (1FAV HHMM-DD). The positions of all sources in the sky are displayed in 
Figure [5} We produced light curves of relative flux variations with FAVA for all sources and made 
them publicly available onlin^j One example light curve is shown for the position of the high-mass 
X-ray binary Cyg X-3 in Figure |6] 

We looked for associations of FAVA sources with previously known variable LAT sources. 
We searched for counterparts within radius R s , which is defined as the 99% statistical error on 
the source position plus the systematic error r sys . R s was deliberately chosen large, to include 
all possible counterparts. In cases for which more than one counterpart is found within R s , we 
consider the closest one. The values of R s for each source and the found counterparts are listed 
in Table [T] We note that the associations were made purely on the basis of positional coincidence. 
We therefore caution that the associated sources should be considered as likely counterparts only. 
For a more confident source association temporal and spectral information need to be considered. 
Additionally, the positional localizations of the FAVA sources could be improved by analyzing each 
source individually with standard likelihood techniques. This is beyond the scope of this paper. 
The associations were assigned as follows: 

• We searched for counterparts among the variable sources in the second Fermi source catalog 



(2FGL, Nolan et al. 2012). We restricted the search to the 458 2FGL sources that have a 
probability of less than 1% of being constant on monthly time scales. We find a variable 
2FGL source within the search radius for 170 of the FAVA sources. For those sources where 
no 2FGL counterpart was found we searched for association with sources in the first Fermi 



source catalog (1FGL, Abdo et al. 2010b). The reason is that sources that flared only once 
at the beginning of the mission might be detected in the 1FGL but not in the 2FGL due 
to the increased integration time in the latter. We restricted the search to the 241 1FGL 
sources that have a probability of less than 1% of being constant, finding an association for 

s https://www-glast.stanford.edu/pub_data/585 



10 



one source. 

• We searched for positional coincidences with Fermi-LAT detected gamma-ray burstsj (GRBs) 
Even though GRBs have typical duration from a few seconds to minutes, their emission is 
sometimes bright enough to be detected over a time scale of 1 week. We find a GRB within 
the search radius for 4 FAVA sources. These FAVA sources flared only once and we have 
verified that the flare occurred during the week of the GRB outburst. 

• We searched for counterparts among LAT sources that were announced in Astronomer's 
Telegrams^ ( ATels) . These sources were found by the automated sky processing (ASP) used 



by the LAT Collaboration (Atwood et al. 2009). We found positional coincidences with 17 



sources. 



We found LAT counterparts for 192 of the 215 FAVA sources. Most of the associated sources, 
177, are AGN. All associations found at higher Galactic latitudes (|6| > 10) belong to this class. 
Among the AGN associations 129 belong to the class of Flat Spectrum Radio Quasars (FSRQs) 
and 29 of them belong to the BL Lacertae (BL Lac) class. The number of gamma-ray emitting 
BL Lacs is approximately the same as the number of FSRQs, FSRQs therefore flare intrinsically 
more frequently in the LAT energy range. This is in agreement with the observation that FSRQs 



are more variable in gamma-rays on monthly time scales (Ackermann et al. 2011b). Three of the 



FAVA sources are associated to non-blazar AGN: two sources are associated to Narrow-Line Seyfert 



galaxies, which were recently found to be variable gamma-ray sources (1FAV J0320+41, Donato &; 



Perkins 2011; 1FAV J0948+01, Foschini et al. 2012), and one source is found coincident with the 



radio galaxy NGC 1275 (1FAV J0320+41, Kataoka et al. 2010). The remaining associated sources 
are AGN of unknown type. 



4. Gamma-ray flares in the Galactic plane 

Of the 215 FAVA sources, 27 are detected at Galactic latitudes smaller than 10°; their positions 
are shown in the lower panels of Figure [5j We found associations to previously known LAT sources 
for all of them. The low-latitude FAVA sources can be grouped into a few categories: 

• Sources associated with Galactic sources: 4 sources coincide with variable 2FGL sources 
which are associated to Galactic sources. These are high-mass X-ray binaries Cyg X-3 (1FAV 
J308+41) and LSI +61 303 (1FAV J0237+60), the Crab Nebula (1FAV J0532+21), and the 
nova V407 Cyg (1FAV J2103+45). In addition, three sources are found coincident with LAT 
sources announced in ATels: the high-mass X-ray Binaries PSRB 1259—63 (1FAV J1303— 64, 



4 http://fermi.gsfc. nasa.gov/ssc/observations/types/grbs/ 



3 https://www-glast.stanford.edu/cgi-bin/pub .rapid 



http://www.asdc.asi.it/feratel (status November 2012) 
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Abdo et al. 2010a 2011a ) and two sources that are likely associated to novas (1FAV J1752— 32, 
Cheung et aL|2012bt 1FAV J0640+05, [Cheung et al.||2012c|d[ ). 



Sources associated with blazars: we find associations with variable 2FGL sources for 15 
sources and to one 1FGL source (1FAV J1124— 64), which were classified either as blazars or 
as AGN of uncertain type in the respective catalogs. Additionally, three sources were found 



coincident with sources announced in ATels and associated to blazars (1FAV J0537+13, Ori 
enti fc D'Ammando||20T2l 1FAV J2202+50, [Ciprini fe Haysl[2012l 1FAV J1718-52, |Chomiuk 



et al. 2013). 



Source with association of unknown type: one source is associated to counterparts of unknown 



type announced in ATels (1FAV J1038— 53, Ciprini et al. 2012). A compact radio sources is 
found coincident with its position. 



Based on these associations 7 FAVA sources are located within the Milky Way. The source 
with the radio counterpart of unknown type might also be Galactic. Additionally, some of the made 
associations with blazars might turn out to be wrong. However, we can already infer statistically 
that most of the sources that we have not associated with Galactic sources are indeed extragalactic 
by calculating the number of expected extragalactic flares at low latitudes. The derivation relies 
on two assumptions: 

1. The majority of the sources at high Galactic latitudes are extragalactic. 

2. Extragalactic sources are isotropically distributed in the sky. 



We derive the number of extragalactic sources within 10° of the Galactic equator from the density 
of sources at latitudes greater than 30° . After considering the difference in solid angle one expects 
41 extragalactic sources at low latitudes. To take into account the reduced sensitivity for flare 
detection in the Galactic plane, we assigned random positions at low latitudes to the high-latitude 
sources. We determined the fraction of flares which would still be detected at the reduced sensitivity 
for each source. This results in an expectation of 24.6 variable extragalactic sources at Galactic 
latitudes smaller than 10°. The probability to detect 20 or fewer flares at low Galactic latitudes is 
21%. The 20 sources not associated with Galactic sources are therefore compatible with being all 
extragalactic. No more than 6 of them can have a Galactic origin at > 90% confidence. 



We note that two gamma-ray binaries LS 5039 and 1FGL J1018. 6-5856 (Abdo et al. 2009a 



Hadasch et al. 2012 ; Ackermann et al. 2012b ) are not detected by FAVA . Their orbital periods of 



3.9 d and 16.6 d, respectively, result in average weekly flux variations below the sensitivity for flare 
detection by FAVA. The X-ray binary Cyg X-l is also not found by our analysis. The flare reported 



from this source by Sabatini et al. (2010) could not be confirmed by the LAT Collaboration^ 



3 http://fermisky.blogspot.de/2010/03/lat-limit-on-cyg-x- l-during-reported.html 
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5. Summary &; Outlook 

We have presented the analysis tool FAVA, which searches for variable gamma-ray emitters 
in Fermi-LAT data. We used FAVA to search for gamma-ray flares on weekly time scales over 
the entire sky. From these flares we derived a list of 215 flaring gamma-ray sources. A list of 
sources, their light curves and their associations are available online^] We searched for positional 
coincidences of these sources with previously known LAT sources, finding counterparts for 192 
sources. We associated 177 sources with AGN and find that FSRQs flare more frequently than BL 
Lacs. 

Twenty-seven of the FAVA sources are located at Galactic latitudes less than 10°. We associ- 
ated 7 of these to known Galactic sources. Among the remaining 20 sources, we found no evidence 
for new gamma-ray transients in our Galaxy. On the contrary, we showed that the majority of 
them are probably extragalactic. For 19 sources we find positional coincidence with AGN. The one 
remaining source is associated to compact radio sources of unknown type. Future multi- waveband 
observations may reveal its nature. 

No flare was detected from a pulsar other than the Crab Nebula. It remains puzzling why this 
source is the only one of its kind to exhibit long-term variability and flaring behavior. We cannot 



confirm the hypothesis reported by Neronov et al. (2012), that gamma-ray variability might be 



common in young pulsars. We also detected no flares associated with previously undetected X-ray 
binary systems. These systems appear to be less efficient gamma-ray emitters than expected before 



the beginning of Fermi observations (Dubus 2007). We note that we detected no flares from the 



Galactic center region, which might have been expected if its gamma-ray emission was linked to 
accretion on the central black hole Sgr A* (Aharonian et al. 2008). 



In the future we plan to apply FAVA on different time scales, and scan the gamma-ray sky for 
short-term flares on time scales of a few hours, and for long-term flux variations of a few months. 
Furthermore, we intend to run the analysis routinely to search for flares as soon as the LAT data 
are processed and sent to the Fermi Science Support Center. This will complement the ASP flare 
search currently used by the LAT Collaboration, and will help to alert the astrophysics community 
about gamma-ray flares in real time. 

Rolf Buhler acknowledges generous support from the Fermi guest investigator program. The 
Fermi LAT Collaboration acknowledges generous ongoing support from a number of agencies and 
institutes that have supported both the development and the operation of the LAT as well as 
scientific data analysis. These include the National Aeronautics and Space Administration and the 
Department of Energy in the United States, the Commissariat a PEnergie Atomique and the Centre 
National de la Recherche Scientifique / Institut National de Physique Nucleaire et de Physique des 
Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in 



7 https://www-glast. stanford.edu/pub_data/585 
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Above 100 MeV 



Above 800 MeV 



CD 



LU 



TO 




Fig. 1. — Illustration of the map generation in FAVA. The week shown covers the time interval 
MJD 54864.655—54871.655 (2009-02-02 to 2009-02-09). The measured counts map is shown in the 
first row, for energies > 100 MeV (left) and energies > 800 MeV (right, white color indicates that 
no counts were detected from this area). The middle panels show the expected counts from the 
average emission observed during the first 47 months of Fermi observations. The generation of 
measured and expected counts is explained in the text. The third row shows the significance of the 
flux variations. As an example the position of the Crab Nebula is indicated by a star in the lower 
left panel. Its flux was increased compared to average with a significance > 5<r for energies > 100 
MeV during this week. The flare is only detected in the low energy range, as the energy spectrum 



was very soft (photon index ~3.5, see Abdo et al. 2011c). An example of a flare detection with a 



negative flux variation is given by the blazar 3C 454.3. Its flux was lower than average for both the 
>100 MeV and >800 MeV energy ranges with a significance > 8a. Figures are shown in Galactic 
coordinates in a Hammer-Aitoff projection. Note that the color scales were adjusted to different 
ranges for the low- and high-energy bands. 
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Fig. 2. — Minimum flux increase required for a flare detection in FAVA in a one- week interval with 
a significance of 5.5a in the low-energy or the high-energy band. Detection above these thresholds 
results in inclusion in the list of flaring sources, as discussed in section [3J The figures are shown 
in Galactic coordinates in a Hammer-Aitoff projection. The flaring source is assumed to have a 
power-law spectrum in energy. The left panel shows the sensitivity for spectrally hard flares (photon 
index 1.5), the right panel for soft flares (photon index 2.5). 
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Fig. 3. — Distribution of the significance of flux variations for the low- (left) and high- energy 
(right) intervals. The integral of both distributions has been normalized to one. Open circles show 
the simulations of a constant gamma-ray sky. Crosses show the measurement of the first 47 months 
of Fermi observations. The solid line in the left panel shows the best-fit Gaussian model for the 
simulated distribution (see text). 
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Fig. 4. — Distribution of angular distances between flares exceeding the significance criteria for 
source detection and known variable gamma-ray sources. The solid and dashed lines show the 
best-fit parametrizations for flares with a detection in the high energy band, and only in the low 
energy band, respectively (see text). 



19 
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1=240 



Fig. 5. — FAVA sources are shown in Galactic coordinates and a Hammer-Aitoff projection. Red 
crosses mark sources for which at least one flare was detected in the high-energy band. Sources that 
were detected only at low energies are marked by a yellow X. The colored image in the background 
shows the maximum significance a m ax detected in each pixel either in the > 100 MeV or > 800 
MeV energy bands during the first 47 months of Fermi observations. The lower three panels show 
the region of Galactic latitude within 10° of the equator, the region enclosed by the dashed lines 
in the upper panel. 
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Fig. 6. — The top panels show the relative variations of counts AN re i = (N — N exp ) / N exp from the 
direction of Cyg X-3 (1FAV J2034+40) in weekly time bins, where N is the number of measured 
counts and N exp is the number of expected counts from the average emission. Bottom panels show 
the significances of these variations. Left-hand plots are for photon energies above 100 MeV, the 
right-hand plots are for photon energies above 800 MeV. 
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